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Abstract 

We present a high statistics study of the leptonic decay constant fp of heavy pseudoscalar 
mesons using propagating heavy Wilson quarks within the quenched approximation, on 
lattices covering sizes from about 0.7 fm to 2 fm. Varying j3 between 5.74 and 6.26 we 
observe a sizeable a dependence of fp when one uses the quark field normalization that 
was suggested by Kronfeld and Mackenzie, compared with the weaker dependence ob- 
served for the standard relativistic norm. The two schemes come into agreement when 
one extrapolates to a ^ 0. The extrapolations needed to reach the continuum quantity 
fp introduce large errors and lead to the value fp = 0.18(5) GeV in the quenched ap- 
proximation. This suggests that much more effort will be needed to obtain an accurate 
lattice prediction for fp. 
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1. Introduction 



The weak decays of D- and 5-mesons will allow to extract the Cabibbo-Kobayashi- 
Maskawa matrix elements from the experimental data, once the hadronic interactions 
in the relevant weak matrix elements are determined from QCD. Lattice calculations in 
principle enable us to compute these QCD effects without model assumptions. 

The simulation of heavy-light quark systems in lattice QCD requires lattice spacings a 
smaller than the inverse relevant masses. At present one reaches values ~ 2 to 4 GeV 
such that D-meson properties are amenable to lattice techniques. 



The 5-meson cannot be handled in such a direct way with lattice methods. The static 
approximation on the lattice 0] enables us to calculate the properties of pseudoscalar 
mesons in the limit of infinite quark mass, and it is very natural to estimate the 5-meson 
matrix elements by interpolating between the Z)-mass range and the static point. The 
latter has been investigated by various groups 0, |^, |[ ^ H, ^ || , and it has been shown 
that (depending on the way one sets the scale) within the range 5.7 < [3 < 6.3, the value 
of the leptonic decay constant fstat changes significantly with the lattice resolution 0, ^ . 
Although the conventional i. e. nonstatic treatment of 'heavy' quarks (applied within — 
and slightly beyond — the charm region) has been pursued by a number of groups | T0| , 
the a dependence of fp has not yet been systematically studied. 



As actual lattice calculations with relativistic heavy quarks do not respect the inequality 
aniq -C 1, lattice artifacts are an important issue. It has been suggested that lattice 
artifacts might be reduced by a modification of the relativistic normalization of states. 



2k, at order 0{amg) p, |TT], |T^, |I3|, |IJ]. In particular, Kronfeld and Mackenzie []T3l have 



given arguments, that heavy quark fields in the Wilson formulation should be normalized 



by yl — 3k/4kc, which differs significantly from the naive normalization. This has been 
shown to yield a much smoother 1/Mp behaviour towards the static point at a fixed lattice 
spacing §, 0. 

The real check of improvement is to observe a flatter a-dependence of fp at fixed rrig. 
In this paper we address this issue by presenting the results of a high statistics study 
of lattices at (3= 5.74, 6.00 and 6.26 (with lattice extensions between 0.7 fm to 2 fm 
and inverse lattice spacings between 1.2 and 3.2 GeV). The results were obtained using 
standard Wilson fermions in the quenched approximation. Smearing techniques were 
needed to improve the saturation of the weak current correlators by the lowest lying 
states at large times. 

The paper is organized as follows: In Section 2 we discuss some details of the calcu- 
lation, in particular the smearing of operators applied to improve the signals at large 
time separations. We present a factorization test among local-local, local-smeared and 
smeared-smeared correlators to ensure ground state dominance. In Section 3 we briefly de- 
scribe the Kronfeld- Mackenzie proposal to normalize heavy-light states. Our main results 
are contained in Section 4 where the a-dependence and finite size effects of the leptonic 
decay constant is presented and the implication of the Kronfeld-Mackenzie normalization 
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at finite values of a is discussed. After extrapolation to the continuum, we perform the 
interpolation between conventional results and the static point, to obtain the estimate for 
fs- Summary and conclusions follow in Section 5. 



2. Operators and Smearing 

The starting point of this paper is our recent study of leptonic decay constants of heavy- 
light mesons in the static quark limit There, we optimized our trial wave functions 
with respect to maximal ground state dominance. In the following we will use the 'best' 
wave function, i.e. a 'gaussian type' wave function with parameters n = 100 and a = 4, 
from ref . 0] , where the reader will find more details about our simulation. The parameters 
of the lattices studied in this work are listed in table 0. 

The pseudoscalar decay constant fp is extracted from the lattice matrix element through 
the relation 

<0\M\%JP>=Z2'./M^fpa'/' . (1) 

Za is the axial current renormalisation constant, ■M.''"^^^ is the fourth component of the 
local axial vector current, and \P > denotes the pseudoscalar ground state. The precise 
normalization is further discussed in section 3. We define a generalized lattice current by 

Mrix,t) = h{x,t)Tl\x,t) (2) 

where 

l'ix,t) = J2^'ix,y,m)mt) (3) 

y 

is a smeared light (/) quark field obtained by applying the trial wave function and 
h{x, t) the local heavy (h) quark field.Q 

The aim is to extract the local matrix element by using wave functions optimized to yield 
early ground state dominance in the meson-meson correlator 

Cl''{t)=Y.<MUx,t)[Mi{0,0)]^> . (4) 

X 

We consider four methods to extract the local matrix element < OlA^^ff^ IP > from the 
local-smeared and smeared-smeared correlators: 

(a) A three-parameter simultaneous fit to Ci°^^f^{t) and C:!^f^^^{t) with equal mass in both 
correlators. 

(b) Same as (a), but with correlators C^'^(t) and Y.x < Mi°^^^{x,t) [A1;^^(0, 0)]^ >. 

Since only smeared light quark fields are needed for the present investigation we suppress the source 
index for the heavy quark field which was used in ref. Q. 
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(c) A constrained fit to the local-smeared correlator Cl^^^^if), using the mass extracted 
from the smeared-smeared correlator, C:!^'^^^^ (t) . 

(d) A fit to the ratio 

(^loc,J (^\ 

R^t) = < ^\M\%^\P > e-''-'/' . (5) 

V C7475 (t) 

All fits (unless stated otherwise) have been carried out with the inverse covariance matrix 
but we have checked that the mean values are not significantly different from the ones 
resulting from a diagonal fit- Ground state dominance has been monitored both by 
and by the plateau in the local mass 



/ip = In — Yj — ~ — cases (a) - (c) (6) 



or 

Bit) 

/^ii(^) = In — r in case (d). (7) 

K(t — a) 

In table we show the results of the different methods for various heavy-light quark 
combinations. Methods (a) - (d) are seen to yield compatible results once ground state 
dominance is established by a low value. 

For the subsequent analysis we chose to apply method (d). The quality of the ground 
state dominance achieved is demonstrated in figure which shows the local masses for 
(3 =5.74, 6.00 and 6.26 at a light quark mass of about 2ms. We hst our results for fp and 
Mp at several combinations of quark masses in table |^. 



2.1 How good are local operators? 



Having established the success of our smearing techniques we are in the position to ex- 
amine the onset of ground state dominance in the case of purely local correlators. The 
ratio of smeared and local correlators [O] 



(Cr''^(t))^ 



rr{t) = ^XL,A.. ^ (8) 



can be expressed in terms of completely local quantities for sufficiently large time sepa- 
rations t, 

( -{m+A)t I -{m+A){T-t) \ 

rrit) = [l + aL I ^_^(,_, J , (9) 

if the correlators are dominated by two states with energies m and m + A: 

Cp' = A? AO(e-"^* + e-'"(^-*)) + A) A;^(e-("^+^)* + e-('"+^)(^-*)) (10) 
= <0|A^;,jP,n> 



and as long as smearing has been operative in the sense that it has suppressed the excited 
state significantly: 

aj ^ A)/A'j « a,oc ^ A]JAI^ . (11) 
The merits of this factorization test are the following: 

• up to terms of order the ratio is purely dependent on the local parameters, 

• the realm of ground state dominance can be clearly identified since the height of the 
plateau is known to be one, 

• the ratio is dimensionless and has a continuum limit in the range where eq. holds. 



In order to determine the onset of ground state dominance in the purely local case we use 
the combined set of the present data and our previous local results [|, Q as an input to 
eq. H We show the ratio ^^475 in figure ^ at /? = 6.26 for fixed light quark mass {mi ~ 2ms) 
and different heavy quark masses as a function of the time separation in units of GeV~^. 
It appears that the plateau regime is reached only at about 5 GeV~^. The local data 
alone seem to allow for an earlier plateau in the local mass. According to the present 
test, however, exploiting the local data under the assumption of ground state dominance 



at smaller time separations leads to an overestimation of fp |[T6|. The pre-asymptotic 
behaviour according to eq. ^ is reflected by the data with a mass gap in the region of 600 
MeV as illustrated by the dashed hues in figure 0. 

At all values of the heavy quark mass that we investigated, the gap is consistent with 
~ 600 MeV. The amplitude aioc does, however, increase significantly when the mass of the 
heavy quark is increased. This entails that the mass- dependence of the decay constant is 
affected by fitting too early. On the basis of the present data, we estimate that this effect 



is of the order of the statistical errors for the (3 = 6.4 data of ref. as well as for the 



P = 6.26 results of ref. |T6 



In the above interpretation we rely on eq. 11. To some extent this equation is checked 



by the lattice spacing independence of rp which results from eq. |TT] (of course, as for any 
observable that has a continuum limit, the normal kind of a-effects should be there). In 
figure we display rp for different values of /?. In agreement with eq. 0, we find that 
is independent of (3 within the statistical errors. 



3. Renormalization 



Our quark fields are normalized according to the relativistic normalization (see ref. ^ for 
the form of the action). Due to the broken chiral symmetry in the Wilson formulation, 
the lattice axial vector current needs a finite renormalization [ITBl 



A^ = ZaM 



loc 

7m 75 



(12) 
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As an intrinsic short distance quantity, Za can be calculated in perturbation theory, which 



at one- loop order gives [19 



Za = 1 - 0.1333/ . (13) 

Since the next order term is so far unknown, it is not a priori clear which coupling con- 
stant g should be used to evaluate Za numerically. It is known that the bare coupling 
go is not a good expansion parameter [^. A recent analysis of a number of short dis- 



tance dominated quantities by Lepage and Mackenzie |20| revealed that a mean field 



approach [^, which absorbs lattice tadpole contributions into effective quantities in the 



Wilson action, leads to a substantially better agreement between lattice Monte Carlo re- 
sults and their perturbative expansions. We therefore use in the following the mean field 
improved coupling 

f = gl/P with P =< ^TrP^, > , (14) 

as determined from the Monte Carlo simulation. In detail, also other forms of mean field 
improvements are being used [^, but we note that they do not differ appreciably in 
numerical values and we have to remember that the errors remain of order 0(^^) in any 
case. 



Kronfeld and Mackenzie have argued that for quark masses that are comparable to 
the inverse lattice spacing, an - at order 0{amq) - different normalization of the axial 
vector current should have matrix elements with smaller lattice artifacts. For comparison, 
we use this nonrelativistic normalization 




with / and h labelling the light and the heavy quark field. Here, Kc denotes the critical 
value of the hopping parameter as determined from the vanishing of the pion mass. The 
renormahzation constant Za is given by [^] 

Za = 1 - 0.0248 f . (16) 

In the following section we will use both and for analyzing the data. The latter is 
denoted by 'KMc norm'. 



4. Results 



4.1 Finite a effects 



In order to determine physical quantities from lattice calculations, it is crucial to perform 
the extrapolation a 0, within a chosen scale. In principle one would like to use f^^ for 
this purpose, since then the 0(^^) uncertainty in Za cancels. At present, the statistical 
errors of /^r in our simulation are however too large to pursue this approach. Therefore we 
will use the string tension, which has small statistical errors, during the extrapolation and 
convert the result finally into the most 'natural' scale, f^^. As stated in the introduction. 
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the supposed improvement of the normahzation ( eq. |T5| ) should manifest itself in a 
smaller slope of /p(a) as a goes to zero [0. The question is whether this is actually the 
case in the available range of f3. In order to check this we have computed fp both with 
the naive (standard relativistic) and the KMc norms at several pseudoscalar masses. The 
data from table ^ are extrapolated in the standard manner to the chiral limit (see ref. [Q) 
and are compiled in table We have determined the a-dependence at fixed values of 
Mp within the range 1.1 GeV < Mp < 2.3 GeV. In figure ^ we show the situation with 
relativistic norm and with the KMc norm. Note, that we had to interpolate between 
the computed Mp values in order to tune for fixed physical Mp as (3 is varied. The 
interpolation introduces only a minor uncertainty since the dependence of fp on Mp is 
very weak. It was estimated by comparing two-point and three-point interpolations and 
added to the statistical errors. 

It can be seen that the KMc norm has a sizeable impact on fp at present values of 
the lattice spacing. Contrary to the expectations, however, fp{a) is turned from a very 
weakly decreasing into a strongly decreasing function. As we increase Mp the effect 
becomes more pronounced. Nevertheless the extrapolated values at a = coincide. We 
have used a linear extrapolation and omitted those points where the effect of the KMc 
factor amounts to more than a 60 % change. The linear extrapolation functions have 
been followed back to the omitted points as dashed lines. These demonstrate that the 
extrapolation is stable, namely, within the statistical errors one would obtain the same 
result at a = if all points were included in the extrapolation. Note that according to 
perturbative arguments, the leading lattice spacing dependence should be linear in a, if 
the 0{g'^) terms are numerically small. We come back to that point in section 4.3. 

In figure ^ we display the effects of the KMc factor onto the quantity 



The first impression is that at fixed and finite values of a the data in KMc norm interpolate 
smoothly between the D region and the static point and this has been interpreted as 
improvement . Given the conclusion from figure ^ the improvement does not survive 
the limit a — 0, however. This is visualized in figure ^ by the error band which refers 
to the results of the above continuum extrapolation. Note also, that the poor scaling 
behavior of fp in the KMc norm is already visible in figure 5.18 of ref. H. 



4.2 Finite size effects 



The data we discussed above was obtained on lattices with periodicity in space L of 
magnitude Ly/a ~ 3. Here we intend to consolidate our results with respect to possible 
finite size effects. For that purpose we have plotted in figure ^ /p as a function of L for 
three pseudoscalars with fixed light quark mass of about 2ms at /5 = 5.74. As expected 
finite size effects decrease with increasing heavy quark mass. 

In the range L^/a > 3 we find no significant finite volume effects with a precision of the 
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order of ~ 2%. Decreasing the light quark mass down to we compare the 8^ X 24 and 
12^ X 24 lattices at (3 = 5.74, c.f. table |. There appears to be a small volume effect 
of maximally 3% ± 2%. This is, however, statistically not significant. So we have kept 
the volume fixed at L^/o^ ~ 3 throughout our final analysis, but estimate from the above 
analysis, that there may be a ~ 5% error involved with this. 



4.3 Conversion to scale /, 



The most 'natural' scale for fp is /vr, since the uncertainty due to the renormalization 
constant Za cancels out in this case and one may also hope, that the ratio fp/fn is 
not affected strongly by the quenched approximation. Lattice measurements of f-j^ are 
generally accompanied by large statistical errors and therefore we have decided to convert 
our results to this scale only after having performed the a extrapolation of fp. To 
achieve this we have decoupled the extrapolations according to 

^(a^O)^ ^^^f'— °| . (18) 

In In/ V(X(a 0) 



To obtain the denominator of eq. |T8| we used both our own data and the results of other 
lattice groups P, H, |2^, |2^ (all data has been changed to the relativistic normalization). 
The compilation of the data is shown in figure |^ Since the a dependence of fj^/ y/a is 
obviously weak, a linear extrapolation to a = is well justified and leads to : 

f 

•'^(a = 0) = 0.269(12) . (19) 
a 



In addition to the quoted statistical error, this quantity has an uncertainty due to the 
missing 0{g^) terms in Za- This should, however, largely cancel out when we take the 
ratio eq. [13. 



4.4 Heavy mass extrapolation 



In figure | we display our final results^ at a = in the form fp{l/Mp), together with our 
static value from ref. 0. As described in the previous section the scale has finally been 
converted to f-,, = 132MeV. The new data appears to depend only weakly on Mp. Note 
however, that the data points carry error bars of order 25% and therefore do not exclude 
a stronger variation in Mp. This is due to the extrapolations to the chiral and continuum 
limits. 



Given this situation we draw an error band that links the conventional results with the 
static point. The Mp dependence of the error band was chosen according to the ansatz 



^To be specific the extrapolation has been performed on the data using the relativistic norm since it 
involves smaller statistical errors than using the KMc norm. 
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At the location of the B meson the error band corresponds to 

fs = 0.18(5)GeV . (21) 

It is evident from figure § that this value is strongly affected by the size and uncertainty 
of /stat- From the direct continuum extrapolation at the mass of the D, we quote 

fn = 0.17(3)GeV . (22) 

At a given Mp one can extrapolate to a — > the ratio fpjfp. In the ratio scale ambiguities 
drop out and additional systematic errors e.g. due to quenching may be reduced. We find 
that the heavy mass dependence of the ratio is very weak and can be taken as constant in 
the mass range between the charm and the B mesons. In the continuum limit we obtain 
fps/fp = 1.09(2) at Ly/a ~ 3, where the statistical error is quite small. We therefore 
have to remember the finite size uncertainty and finally quote 

/p7/p = 1.09(2)(5). 



5. Summary and conclusions 

We have carried out a high statistics analysis of the pseudoscalar decay constant. In the 
case of purely local operators, the ground state signal is attained only at time separations 
larger than 5 GeV~^. For this reason smearing techniques should also be applied in the 
nonstatic situation. 

We observe little volume dependence of fp as long as L^/a > 3. For 5.74 < P < 6.26, we 
discover a marked variation with a, if the KMc normalization is used, while the relativistic 
normalization results in /p-values that are rather insensitive to our variation of a. This 
implies that KMc does not yet lead to improvement in the sense that finite a results are 
shifted towards their continuum limit. 

In this context, we emphasize that it is misleading to look only at the Mp dependence 
of fp at one value of jS. One rather has to work at different f3 values and to extrapolate 
to the continuum. Proceeding in this manner we find that the difference originating from 
the two normalizations vanishes, as it should. 

We consider the present investigation as exploratory in spite of the statistics of about 
100 independent configurations on the largest lattice. The estimates of fp and Jd (see 
eqs. ^,^) still carry rather large errors originating from the unavoidable extrapolation 
in a. Study of the a dependence is also necessary when an improved action is used. The 
extrapolation errors can be greatly reduced with the power of available parallel super- 
computers, by increasing the number of points in a with very high statistics. 

As we have mentioned, it would also be desirable to perform the whole continuum ex- 
trapolation on the ratio fp/ f-n in order to be completely free of the uncertainties in Za- 
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It appears, however, that better ways of computing /,r need to be found before this is 
feasible. 

Along the same lines, a study of the finite a effects of Bb is necessary, since it has been 
clearly demonstrated for the case of that a effects do not necessarily cancel out for 



these observables [24 
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Tables 



Table 1: 

Pseudoscalar decay constant and mass in lattice units as a function of tmin, the small- 
est time separation used in the fit. Smearing is used for the light quark propagator. 
The results correspond to a light quark mass of about twice the strange quark mass (or 
M{l,lf/a ~ 4). denotes xVd.o.f.. 







P = 5.74 Ki = 0.1560 


K2 = 
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1.353(4) 


2.9 






3 
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1.080(2) 


2.0 






4 


0.1277(14) 


1.079(2) 


2.0 






5 


0.1260(15) 


1.078(3) 


1.0 






6 


0.1266(15) 


1.078(3) 


1.1 






7 


0.1273(19) 


1.079(3) 


1.1 






9 


0.1270(19) 


1.079(4) 


1.4 






18^ X 36 9 conf. 






^min 


af/ZA 


aM 


x' 






6 


0.1281(54) 


1.084(7) 


uncorr. 






7 


0.1264(57) 


1.082(7) 


uncorr. 






9 


0.1208(76) 


1.075(11) 


uncorr. 





f3 = 6.26 Ki = 0.1492 K2 = 0.1200 





12^ X 48 


103 conf. 




^min 


af/ZA 


aM 


x' 


4 


0.0859(23) 


1.040(5) 


1.8 


6 


0.0845(34) 


1.039(8) 


2.0 


8 


0.0834(30) 


1.039(6) 


1.6 


10 


0.0853(46) 


1.039(7) 


1.2 


12 


0.0831(47) 


1.035(7) 


1.2 


14 


0.0808(48) 


1.031(9) 


1.2 


16 


0.0823(55) 


1.033(8) 


1.5 




18^ X 48 


76 conf. 




^min 


af/ZA 


aM 


x' 


4 


0.0833(21) 


1.039(7) 


2.9 


6 


0.0813(21) 


1.034(6) 


2.0 


8 


0.0811(27) 


1.034(5) 


2.3 


10 


0.0784(29) 


1.030(5) 


1.3 


12 


0.0777(29) 


1.029(5) 


1.5 


14 


0.0805(39) 


1.033(6) 


1.5 


16 


0.0783(44) 


1.030(7) 


1.7 



Table 2: 



af /Za calculated with different methods as described in the text: 

(a) simultaneous fit to C'^ii and Ci,°!;'/, 

(b) simultaneous fit to and < ^^^^^..l^, ^) [^^4,(0, O)]^ >, 

(c) determine amp by fit to C:^f^^ and use this as constraint in the fit to C*^"^^, 

(d) ratio method. 

We show the results for the lattices used to obtain our final results, actually denotes 
xVd.o.f.. 

/3 = 5.74 Ar| X ATy = 8^ X 24 = 0.1560 ^2 = 0.1250 



^min 


(a) 


(b) X' 


(c) X' 


(d) X' 


2 
3 

4 
5 
6 
7 
8 


0.2274(14) 2.8 
0.2256(20) 2.6 
0.2224(24) 2.1 
0.2203(27) 1.9 
0.2159(35) 1.7 
0.2145(41) 2.0 
0.2181(44) 1.7 


0.2372(26) 8.8 
0.2313(28) 6.9 
0.2250(28) 4.1 
0.2216(29) 3.4 
0.2158(31) 1.4 
0.2142(35) 1.4 
0.2130(38) 1.8 


0.2314(17) 6.8 
0.2280(23) 3.0 
0.2238(28) 2.2 
0.2196(30) 1.6 
0.2169(33) 2.0 
0.2154(36) 2.6 
0.2144(38) 1.3 


0.2270(19) 3.1 
0.2250(19) 3.1 
0.2222(22) 2.5 
0.2207(28) 2.7 
0.2169(35) 2.6 
0.2157(38) 2.5 
0.2209(43) 2.3 


(5 = 6.00 NIxNt^ 12^ X 36 = 0.1525 ^2 = 0.1250 


^min 


(a) x' 


(b) X' 


(c) X' 


(d) X' 


3 
4 
5 
6 
7 
9 
11 


0.1267(14) 1.8 
0.1266(10) 1.9 
0.1251(14) 1.2 
0.1257(17) 1.1 
0.1262(22) 1.2 
0.1248(20) 1.3 
0.1236(26) 1.4 


0.1333(14) 4.9 
0.1318(15) 3.5 
0.1307(15) 1.7 
0.1296(18) 1.7 
0.1282(20) 1.5 
0.1257(23) 1.1 
0.1243(23) 1.1 


0.1309(16) 1.8 
0.1303(16) 1.8 
0.1268(20) 1.4 
0.1281(16) 1.3 
0.1286(20) 1.5 
0.1257(26) 1.8 
0.1231(30) 2.0 


0.1280(17) 2.0 
0.1277(14) 2.0 
0.1260(15) 1.0 
0.1266(15) 1.1 
0.1273(19) 1.1 
0.1270(19) 1.4 
0.1258(24) 1.7 


/3 = 6.26 NIxNt^ 18^ x 48 = 0.1492 ^2 = 0.1200 


^min 


(a) X' 


(b) x' 


(c) X' 


(d) x' 


4 

6 

8 

10 

12 

14 

16 


0.0819(30) 2.5 
0.0809(33) 2.6 
0.0799(30) 2.6 
0.0770(25) 1.6 
0.0775(30) 1.6 
0.0800(36) 1.5 
0.0789(35) 1.5 


0.0842(42) 3.3 
0.0825(39) 2.7 
0.0808(32) 2.8 
0.0782(22) 1.3 
0.0783(24) 1.2 
0.0795(29) 1.1 
0.0793(29) 1.2 


0.0821(36) 2.2 
0.0810(36) 2.3 
0.0811(37) 2.4 
0.0777(33) 1.4 
0.0782(38) 1.7 
0.0800(34) 1.8 
0.0782(41) 1.7 


0.0833(21) 2.9 
0.0813(21) 2.0 
0.0811(27) 2.3 
0.0784(29) 1.3 
0.0777(29) 1.5 
0.0805(39) 1.5 
0.0783(44) 1.7 



Table 3: 



Pseudoscalar decay constant and mass in lattice units. The data has been used in the 
detailed analysis to extract fs and f^, with the exception of the results at /5 = 5.74, 
12^ X 24. Smearing has been applied to the light quark propagator. Ki corresponds to the 
light quark, actually denotes x^/d.o.f.. 



(3 = 5.74 Nix Nt = 8^x 24 





1^2 


fit range 


aM 






0.156 


0.09 


5-12 


1.999(4) 


0.1905(33) 


1.0 




0.125 


4-12 


1.344(4) 


0.2206(24) 


2.8 




0.140 


4-12 


1.065(4) 


0.2231(31) 


3.0 




0.150 


4-12 


0.870(4) 


0.2164(41) 


2.4 




0.156 


4-12 


0.748(5) 


0.2008(57) 


2.4 


0.162 


0.09 


5-12 


1.924(8) 


0.1750(40) 


1.0 




0.125 


4-12 


1.261(2) 


0.2022(23) 


1.8 




0.140 


4-12 


0.969(5) 


0.2037(33) 


2.1 




0.150 


4-12 


0.761(4) 


0.1972(48) 


2.0 




0.162 


5-12 


0.474(8) 


0.1608(58) 


1.1 


0.1635 


0.09 


5-12 


1.904(10) 


0.1697(43) 


0.8 




0.125 


4-12 


1.242(7) 


0.1973(24) 


2.0 




0.140 


4-12 


0.948(5) 


0.1997(30) 


1.1 




0.150 


4-12 


0.734(4) 


0.1917(44) 


1.3 




0.1635 


5-12 


0.392(12) 


0.1479(64) 


1.1 




P = 


5.74 Nl 


xNt = 12^ 


X 24 






1^2 


fit range 


aM 


a I/Za 


x' 


0.156 


0.09 


5-12 


2.001(6) 


0.1872(24) 


0.6 




0.125 


4-12 


1.354(5) 


0.2179(24) 


0.8 




0.140 


4-12 


1.073(7) 


0.2207(12) 


1.1 




0.150 


4-12 


0.878(5) 


0.2133(23) 


1.1 




0.156 


4-12 


0.755(5) 


0.1959(24) 


0.8 


0.162 


0.09 


5-12 


1.927(17) 


0.1689(40) 


0.8 




0.125 


4-12 


1.267(7) 


0.1974(29) 


0.6 




0.140 


4-12 


0.976(9) 


0.2001(24) 


1.0 




0.150 


4-12 


0.768(4) 


0.1929(22) 


1.1 




0.162 


5-12 


0.481(2) 


0.1628(30) 


0.8 


0.1635 


0.09 


5-12 


1.911(15) 


0.1631(47) 


0.7 




0.125 


4-12 


1.246(5) 


0.1914(27) 


0.5 




0.140 


4-12 


0.953(8) 


0.1935(23) 


0.9 




0.150 


4-12 


0.740(2) 


0.1864(30) 


1.1 




0.1635 


5-12 


0.395(4) 


0.1523(35) 


0.9 



(3 = 6.00 Nix Nt = 12^ X 36 





«2 


fit range 


aM 






0.1525 


0.10 


5-13 


1.574(7) 


0.1060(15) 


0.7 




0.115 


5-13 


1.279(3) 


0.1186(14) 


0.6 




0.125 


5-18 


1.077(2) 


0.1260(15) 


1.0 




0.135 


5-18 


0.870(3) 


0.1349(13) 


0.6 




0.145 


7-18 


0.649(1) 


0.1301(23) 


0.7 




0.1525 


4-11 


0.452(2) 


0.1090(15) 


0.6 


0.1540 


0.10 


5-13 


1.546(8) 


0.1003(15) 


0.7 




0.115 


5-13 


1.252(5) 


0.1123(15) 


0.6 




0.125 


5-18 


1.050(2) 


0.1190(13) 


1.0 




0.135 


5-18 


0.840(4) 


0.1243(14) 


0.6 




0.145 


5-18 


0.617(1) 


0.1211(18) 


0.9 




0.154 


4-11 


0.368(3) 


0.0980(16) 


0.7 


0.1558 


0.10 


5-13 


1.521(8) 


0.0925(17) 


1.1 




0.115 


5-13 


1.221(5) 


0.1036(16) 


1.1 




0.125 


5-18 


1.018(5) 


0.1100(17) 


1.2 




0.135 


5-18 


0.806(4) 


0.1149(18) 


0.9 




0.145 


5-18 


0.579(4) 


0.1111(21) 


0.9 




0.1558 


1 11 


0.257(8) 


0.0821(22) 






/5 = 


6.26 N^x Nt = 18^ 


X 48 








fit range 


aM 


a f/ZA 


x' 


0.1492 


0.09 


10-18 


1.639(8) 


0.0577(26) 


0.5 




0.10 


10-18 


1.433(7) 


0.0643(25) 


0.4 




0.120 


10-23 


1.033(5) 


0.0784(29) 


1.3 




0.135 


10-23 


0.707(4) 


0.0880(25) 


1.4 




0.145 


6-23 


0.461(3) 


0.0882(19) 


1.3 




0.1492 


5-18 


0.344(4) 


0.0759(33) 


1.0 


0.1506 


0.09 


10-18 


1.614(11) 


0.0517(28) 


0.9 




0.10 


10-18 


1.409(9) 


0.0578(26) 


0.7 




0.120 


10-23 


1.003(7) 


0.0709(29) 


1.4 




0.135 


10-23 


0.676(4) 


0.0799(24) 


1.3 




0.145 


6-23 


0.424(4) 


0.0810(21) 


1.2 




0.1492 


4-18 


0.301(4) 


0.0672(37) 


1.1 




0.1506 


4-18 


0.254(5) 


0.0642(31) 


1.0 


0.1514 


0.09 


10-18 


1.596(14) 


0.0478(29) 


0.8 




0.10 


10-18 


1.391(11) 


0.0532(28) 


0.7 




0.120 


10-23 


0.985(9) 


0.0662(30) 


1.3 




0.135 


10-23 


0.658(4) 


0.0745(26) 


1.3 




0.145 


6-18 


0.407(4) 


0.0761(21) 


1.1 




0.1492 


4-18 


0.275(4) 


0.0638(39) 


1.2 




0.1506 


4-17 


0.224(5) 


0.0596(41) 


0.9 




0.1514 


4-17 


0.189(5) 


0.0565(39) 


0.7 



Table 4: 



We give the pseudoscalar decay constant and mass in lattice units extrapolated to and 
Ks, as well as the ratio of the decay constant for a light quark fixed to the strange quark 
mass to that evaluated at the chiral limit, for all the heavy quarks Hh considered. The 
subscript u and s denote quantities evaluated at the chiral limit and at the strange quark 
mass respectively. The strange quark mass was fixed using the a scale. 





p — 0.(4 


= 0.1664(5) = 


n 1 cifiOf'Q^ 
U. iooyl^o ) 










a fs/ZA 


a IVlg 


f / f 

J SI J u 


U.UO 


U.iiy ( (^iUzj 


2.502(13) 


0.1313(74) 


z.54z(^y j 


i.Uy / [66 j 


u.uy 


0.iD29(^o2 ) 


1.871(13) 


0.1726(39) 


1.916(9 ) 


1.0oy(^ll ) 


U.izo 


O.io90 (^OO J 


1.205(7) 


0.2001(32) 


i.Zoo[Z ) 


l.Ooo^iOj 


n 1 /I n 
U.i4U 


U. iyU i \ S(i ) 


0.904(6) 


0.2020(29) 


U.yoU(^4 ) 


1 nc;n('i n^ 


0.150 


0.1829(53) 


0.684(5) 


0.1947(42) 


0.749(4) 


1.065(10) 




Q a nn 


Kc = 0.1572(4) Ks = 


U.i£)4U(^o ) 










a fs/ZA 


a IVlg 


f 1 f 

J SI Ju 


O.iO 


\).\)olo\V I ) 


1.498(10) 


0.0936(16) 


1.02o(9 ) 


1.0/0(^1 ) 


O.iio 


0.09oo(^lo ) 


1.197(7) 


0.1048(16) 


i.zzo(4j 


1.0/o(^0J 


O.izo 


O.lOool^lo j 


0.995(5) 


0.1113(16) 


l.Ozoi^z j 


1.0M(^6j 


O.ioo 


0.10oo(^19 j 


0.780(7) 


0.1163(16) 


O.olOl^o j 


1.0 /Z[6 ) 


n 1 /1 1^ 
U.i40 


n 1 nQO^'Qi ^ 

U. l\JOZ,\OL ) 


0.551(2) 


0.1123(20) 


n ^QA( A\ 
U. 084(^4 ) 


1 nss ('1 A \ 
i.Uoo(^i4 j 




(3 = 6.26 


Kc = 0.1524(6) Ks = 


0.1507(4) 




K-h 


a U/Za 


aMu 


a fs/ZA 


a Ms 


fs/fu 


0.09 


0.0437(32) 


1.579(19) 


0.0474(30) 


1.595(13) 


1.086(18) 


0.10 


0.0486(31) 


1.375(14) 


0.0529(29) 


1.390(11) 


1.087(15) 


0.120 


0.0609(33) 


0.965(10) 


0.0656(30) 


0.983(8) 


1.077(13) 


0.135 


0.0689(29) 


0.636(6) 


0.0741(25) 


0.655(4) 


1.074(12) 


0.145 


0.0711(24) 


0.382(4) 


0.0757(21) 


0.403(3) 


1.065(10) 


0.1492 


0.0580(50) 


0.245(4) 


0.0628(42) 


0.272(3) 


1.082(26) 



Table 5: 

The values of the inverse lattice spacing from the string tension using ^/a^ = A20MeV. 
For comparison we also give the values determined from our and rrip data. 



p 






"-nip 


5.74 
6.00 
6.26 


1.118(9) 
1.876(19) 
2.775(18) 


1.35(8) 
2.36(8) 
3.46(37) 


1.41(4) 
2.15(10) 
2.94(16) 



Figure Captions 



1. Local masses as defined in eq. |^ for tlie pseudoscalar meson are plotted versus the 
time separation t in lattice units. Shown are the data for (3 =5.74, 6.00 and 6.26 for 
a light quark mass of about 2'ms and a heavy quark mass corresponding to Kh- 

2. Ratio of smeared and local correlators as defined in eq. ^ at /? = 6.26 for different 
heavy quark masses and a light quark mass of about 2ms. Deviations from unity 
signal contamination of local correlators by excited states. The scale has been set 
by use of f^,, cf. table ||. 

3. Ratio of smeared and local correlators as defined in eq. ^for f3 = 6.0 and 6.26 and 
a heavy quark mass of about the charm quark mass. 

4. fp vs a in physical units for different meson masses fixed by interpolating between 
the calculated ones. The lines are linear fits to the data points using only the values 
where the ratio KMc norm to standard norm is smaller than 1.6. The scale is set 
by the string tension. 

5. fp is shown as a function of 1/Mp at j3 =5.74, 6.00 and 6.26, where the light 
quark mass is extrapolated to the chiral limit. The set of data points labeled by 



"KMc" were obtained with the normalization eq. |15| and the dashed lines denote 
the error band after taking the continuum limit of these data. The others denoted 
by "Standard" are with the relativistic normalization. The static points are taken 
from ref. [Q. The scale was set by y/a. 

6. jpl\fo is shown vs La/o" at /5 = 5.74 for different heavy quark masses. The light 
quark mass was fixed at about 2ms. 

7. Compilation j^j^ as function of The string tension values are as listed 
in ref. with linear interpolation in ln(a^(T). The symbols refer to the results of 
several groups: 

Filled circle: refs. [^, Filled square: this work; Filled triangle: ref. [Q; Filled 



star: refs. [§, H and this work; Open circle: ref. ^Sf; Open square: ref. ||2^; Open 



triangle: this work; Open rhomb: ref. [^]; The star refers to the linearly extrapolated 
value. 

8. /p as function 1/Mp. The scale is taken from The lines correspond to the error 
band described in the text. 
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